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Abstract: Scattering processes in the cosmic microwave background limit the propagation of 
ultra high energy charged particles in our Universe. For extragalactic proton sources resonant 
photopion production results in the famous Greisen-Zatsepin-Kuzmin (GZK) cutoff at about 
4 X 10^® GeV expected in the spectrum observed on Earth. The faint flux of ultra high energy 
cosmic rays of less than one event per year and cubic kilometer and the large systematic 
uncertainties in the energy calibration of cosmic ray showers is a challenge for cosmic ray 
observatories and so far the GZK cutoff has not been unambiguously confirmed. We have 
investigated the possibility that the primaries of super-GZK events are strongly interacting 
neutrinos which are not subject to the GZK cutoff. For the flux of protons and neutrinos from 
extragalactic optically thin sources and a flexible parameterization of the neutrino-nucleon 
cross section we have analyzed the cosmic ray spectra observed at AGASA and HiRes taking 
also into account results from horizontal events at AGASA and contained events at RIGE. 
We find that scenarios of strongly interacting neutrinos are still compatible with the data 
requiring a steep increase of the inelastic neutrino-nucleon cross section by four order of 
magnitude within one energy decade compared to the Standard Model predictions. We also 
discuss the impact of the preliminary cosmic ray spectrum observed by the Pierre Auger 
Observatory. 
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1. Introduction 


The origin and chemical composition of cosmic rays with energies above 10® GeV - so called 
ultra high energy (UHE) cosmic rays (CRs) - is still an open question in astroparticle physics. 
The large-scale isotropy |j^, ^ and also the shower characteristics Q of these events give 
reason to believe that CRs around the ankle at about 10^® GeV are dominated by extra- 
galactic protons. As charged particles, these protons are subject to scattering processes in 
the intergalactic photon gas. It has been shown Q (see also Refs. H) 1) 0) H) tha-t the CR 
spectrum between 10® GeV and 4 x 10^® GeV is in a remarkably good agreement with the 
propagated flux from extragalactic proton sources using a simple injection spectrum with two 
free parameters. In this case, the ankle can be identihed as an e'’'e“ pair production dip of 
protons scattering off photons of the cosmic microwave background (CMB) together with a 
pile-up of protons due to photo-pion production. 


Soon after the discovery of the CMB in 1965 it was pointed out by Greisen Zatsepin, 
and Kuzmin |10| that resonant photopion production of protons above 4 x 10^® GeV limits 
their propagation to about 50 Mpc.^ For an extra-galactic proton source this results in a 
cut-off expected in the CR spectrum. So far the experimental resolution for this feature is 
very poor due to low statistics and large systematic errors in energy calibration (cf. Fig. ||). 
The next generation cubic-kilometer-sized detectors, notably the Pierre Auger Observatory 
(PAO) will soon have the necessary exposure to resolve this issue. 


A continuation of the CR spectrum beyond the GZK cutoff to extremely high energies 
(EHFs) seems only consistent with a proton dominance if the sources he in our local cosmic 
environment. Only very few astrophysical accelerators can achieve the necessary energies (see 
e.g. Ref. [11| for a review) and so far none of the candidate sources have been confirmed. It 
has been speculated that decaying superheavy particles, which could be some new form of 
dark matter or remnants of topological defects, could be a source of UHE CRs, but also these 
proposals are not fully consistent with the gamma ray spectrum at lower energies [p!^ ]. 

These missing parts of the CR puzzle has led to speculations about a different origin 
of EHE CRs. In particular, a flux of neutral components would not be attenuated by the 
CMB and could be fueled by very distant sources. In the late 60s Berezinsky and Zatsepin |]T] 
proposed that cosmogenic neutrinos produced in photopion production of protons in the CMB 
could explain EHE events assuming a strong neutrino-nucleon interaction. The realization 
of such a behavior has been proposed abundantly in scenarios beyond the (perturbative) 
Standard Model (SM): e.g. arising through compositeness [14, 15, 16|, through electroweak 
sphalerons jU, |T^, |l^, |2^, H, ^ |5|, ^j, through string excitations in theories with a low 
string and unihcation scale [25, 26, 27], through Kaluza-Klein modes from compactified extra 

^], or through p-brane production in models with warped 


dimensions 


extra dimensions [34 


I, respectively (for recent reviews, see Refs. [37, 


^For heavier nuclei photo-disintegration in the CMB give a comparable or even stronger attenuation above 
this energy. 
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Figure 1: Cosmic ray spectrum observed by various experiments. Also shown is the size of the 
systematic error in energy calibration. 


Up to now, UHE cosmic neutrinos have been searched for in the Earth atmosphere 
(Fly’s Eye Q and AGASA Q), in the Greenland (FORTE Q) and Antarctic ice sheet 
(AMANDA and RIGE |^]), in the sea/lake (notably BAIKAL |^), or in the regolith of 
the moon (GLUE [^). A large neutrino-nucleon cross section is not necessarily in conflict 
with the very rare UHE neutrino events claimed by these experiments, if the transition from 
the weak SM interaction to a strong contribution of new physics proceeds over a sufficiently 
short energy range. Those neutrinos interacting very strongly will then be hidden in the 
nucleonic CR background . In our quantitative analysis of scenarios with strong neutrino- 
nucleon interactions we have used the search results from AGASA and RICE. 

In this paper we will report on a recent statistical analysis ||4^ investigating scenarios with 
strong neutrino-nucleon interactions as an explanation of GZK excesses. For the flux of CRs 
and neutrinos from distant optically thin sources introduced in §2 and a model for a strong 
neutrino-nucleon interaction introduced in §3 we will show in §4 the results of goodness-of-fit 
test for existing CR data from AGASA |^] and HiRes [^, M] as well as preliminary 

data from PAO |^, We also include the search results on neutrino-induced events from 
AGASA and RICE. Finally, we conclude in §5. 


2. Extragalactic Particle Production at Ultra High Energy 

The hypothesis of extragalactic proton dominance around the ankle is in very good qualitative 
agreement with the observed CR spectrum between 10® GeV and 4 x 10^® GeV I, & i g i 
as shown in Fig. One can distinguish two mechanisms for the production of these UHE 
protons. They could be produced in the decay of superheavy particles, formed in the decay 
of topological defects or being a part of the dark matter content in our Universe. In general. 
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Figure 2: The flux of extragalactic protons (solid line) from Ref. 0 using a power-like injection 
spectrum with red-shift evolution of the source luminosity oc (1 -I- z)'^E~^. Also shown are the fluxes 
of cosmogenic neutrinos (dashed line) summed over all flavors. 


these top-down scenarios involve physics beyond the SM like Grand Unified Theories. In the 
more conservative bottom-up approach charged particles are magnetically confined in large 
astrophysical objects, such as active galactic nuclei or gamma ray bursts and accelerated 
through (repeated) scattering off plasma shock fronts. Photopion interactions with the photon 
gas inside the accelerator produces high energy neutrons, neutrinos, and gamma rays which 
can escape the magnetic region (see Fig. We will assume in the following that extragalactic 
UHE protons are produced in optically thin astrophysical accelerators, which allow for a simple 
determination of neutrino fluxes relative to the flux of CRs. 

2.1 Extragalactic Protons and Cosmogenic Nentrinos 

For an efficient acceleration of a charged particle by repeated scattering processes within 
an astrophysical source it is necessary that its gyroradius is smaller than the size of the 
accelerator. This simple geometric argument is conveniently depicted as the “Hillas diagram” 
shown in Fig. In the case of non-cyclic one-shot acceleration scenarios a similar expression 
for the maximal energy holds From Fig. Q one can see that only very few objects are 

able to produce UHE protons and even less of them EHE protons as candidates for super- 
GZK events. In our statistical analysis we have assumed that UHE protons are injected 
by spatially homogeneous and isotropic distributed astrophysical accelerators with red-shift 
dependent luminosity according to 

^cr{z,Eck) cc {1-\-z)'^ (2.1) 

^min Z <C. 2^max 
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We exclude nearby (redshift Zmin) and early (zmax) sources and fix these parameters in the 
following at Z m\r, = 0.012 (corresponding to r min ~ 50Mpc) and Z max = 2. The maximal 
injection energy illmax is fixed at 10^^ GeV in our analysis which is compatible with the 
highest energy proton sources shown in the Hillas diagram Fig. 

The flux of UHE protons from distant sources is subject to energy redshift, e"’“e“ pair 
production, and photopion-production in the CMB (see e.g. [^, 56, which can be taken 
into account by means of propagation functions The resonantly produced pions, 

which are also the reason for the GZK cut-off, decay into electron and muon neutrinos via 
the reaction chain p + 7cMB —> A^ + 7r+ —> ... —> e~^ ... and provide a guaranteed 

source of cosmogenic UHE neutrinos observed at Earth (see Eig. |^. The initial ratio of 
neutrinos '■ '■ Vt as 1:2:0 is then redistributed as 1:1:1 by oscillation effects. 

2.2 Neutrinos from Optically Thin Sonrces 

The inelastic scattering of the beam protons off the ambient photon gas in the source will 
also produce photopions which provide an additional source of UHE neutrinos (see Eig. 

In general, the relative production of other (neutral) particles in the source depend on details 


such as the densities of the target photons and the ambient gas |^|. In the case of optically 
thin proton sources we can estimate the flux of UHE neutrinos as 


Cy{z,Ey) 


3/:^)£cr 


\ E, 


^ / 




E, 


( 2 . 2 ) 


For the mean energy ratio we have used (Ecr/E^) ~ 0.07, suitable for resonant photopro¬ 
duction at the energies in question (see Ref. ||^]). This flux of neutrinos is directly associated 
with the injected flux of GRs and serves as a minimal contribution from the source. In 




Figure 3: Left panel: Sketch of an optically thin CR source. Each UHE proton is accompanied by 
three UHE neutrinos. Right panel: In optically thick CR sources neutrons may scatter of the photon 
gas before escaping the source. This produces a larger neutrino to CR ratio. 
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Figure 4: Hillas diagram (from Ref. |Q) for objects with characteristic radius R and magentic field 
B. The dashed-dotted lines correspond to different proton energies from 10® GeV to 10^® GeV. The 
solid lines correspond to limits of the maximal proton energy including energy losses and interactions 
for the models discussed in Ref. [|^. The sources are: neutron stars (ns), white dwarfs (wd), sunspots 
(ss), magnetic stars (ms), active galactic nuclei (ag), interstellar space (is), supernova remnants (sn), 
galactic disk (d), halo (h), radio galaxy lobes (rg), clusters of galaxies (cl) and intergalactic medium 
(ig). Also shown are jet-frame parameters for blazars (bl) and gamma ray bursts (gb). 


optically thicker sources neutrons may undergo photohadronic interactions before escaping 
the confinement region as it is sketched in the right panel of Fig. This will decrease the 
emissivity of neutrons compared to that of neutrinos. Depending on the ambient gas, pp 
interactions may also dominate over photohadronic processes in the source and produce ad¬ 
ditional neutrinos. As it was discussed in Ref. 0 this flux of high energy neutrinos from 
optically thin sources is almost in reach of the AMANDA-II detector for neutrino energies of 
the order 10^ GeV and should soon be observable by its successor experiment IceCube (see 
Fig. ^). This observation is crucial for our assumption, that extragalactic protons start to 
dominate in the CR spectrum already at a low crossover below the ankle. If even cosmogenic 
neutrinos are not detected in future experiments, also the hypothesis of extragalactic proton 
dominance above the ankle has to be reconsidered. 


3. Strongly Interacting Cosmic Neutrinos 

The SM interaction length of neutrinos exceeds the vertical column depth by several orders 
of magnitude. For a reasonable contribution of extragalactic neutrinos in vertical CRs one 
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Figure 5: The neutrinos from optically thin sources (tilted solid lines) are close to be measured by 
neutrino telescopes such as AMANDA/IceCube (from Ref. Q). The fluxes are normalized to the 
AGASA and HiRes cosmic ray spectrum, respectively. Also shown are the Waxman-Bahcall flux 
(horizontal solid line) assuming a CR injection as E~^ and the fluxes of cosmogenic neutrinos (dashed- 
dotted lines). 


has to assume a strong enhancement of the neutrino-nucleon interaction. However, the tran¬ 
sition from a weak to a strong interaction has to be sufficiently fast in order to reproduce 
the low neutrino rates or negative search results from experiments such as AGASA ||40t|, 


AMANDA 1^, ANITA BAIKAL Fly’s Eye i3 
RICE [^]. Further ambitious projects include IceCube 
by EUSO H], OWL ||, or SalSA ||]. 


, EORTE GLUE and 
and PAO |^|, possibly followed 


In the introduction we have already given examples of scenarios beyond the (perturbative) 
SM which predicts a strong neutrino-nucleon interaction. In the following we will use a 
flexible parameterization of a strong neutrino-nucleon inelastic cross section (u))^) focusing 
on three characteristic parameters: (i) the energy scale Eth of the new underlying physics, (ii) 
the amplification A compared to the SM predictions, and (iii) the width B of the transition 
between weak and strong interaction. A mathematical convenient parameterization illustrated 
in Eig. ^ is given by 


logio 


a 


uN 


A 


a' 


SM 

vN 


1 

2 


logg 

Eth 


(3.1) 


In general, experiments distinguish different CR primaries by their characteristic shower de¬ 
velopment in matter. Eor the sake of simplicity we will assume that the characteristics of the 
showers induced by strongly interacting neutrinos are indistinguishable from those induced 
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by protons. In particular, we assume for both primaries (i) a complete conversion of the 
incident energy into the shower, and (ii) equal detection efficiencies at the highest energies. 
Under these assumptions the number of detected events can be expressed as 

IVobs = J dES{E){Jp{E)+n{a^j:^{E))ME)) , (3.2) 

where we define as the relative experimental exposure to strongly interacting neu¬ 

trinos compared to the exposure £{E) to protons. In general, the value of TZ{cr'^ff (E)) is 
determined by the search criterion on the zenith angle 6 and the (observed) atmospheric 
depth adopted by each experiments. Figure ^ gives the relative exposure for vertical showers 
at AG AS A and HiRes (see Ref. [^] for details). 


4. Quantitative Analysis 


For the vertical CR spectrum of AGASA, HiRes-I/II, and the Fly’s Eye Stereo data together 
with the search results of horizontal events at AGASA and contained events at RIGE we 
show in the left panel of Eig. ^ the 90%, 95%, and 99% confidence level (GL) of the inelastic 


cross section (Eq. 3.1) from a goodness-of-fit test using the fluxes from optically thin sources 
(Eq. and ^.21) . The details of the statistical analysis and the approximations involved can 
be found in Ref. [46|. The goodness-of-fit test of our model using the combined data requires, 
to the 90% GL, a steep increase by an amplihcation factor of A > 10^ over a tiny energy 
interval B < 10^'® at about 10^^ GeV. Within this particular scenario of strongly interacting 
neutrinos, the SM neutrino-nucleon inelastic cross section corresponding to A = 1 is not 
favored by the data. 




E [GeV] [mb] 


Figure 6: Left panel: An example of the strong neutrino-nucleon interaction given in Eq. &■ 
Right panel: Relative exposure TZ (Eq. (^)) to strongly interacting neutrinos as a function of the 
neutrino-nucleon inelastic cross section of AGASA (solid) and HiRes (dashed), respectively). 
The hatched region shows the predicted contribution from SM charged current interaction [M M. 
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Figure 7: The range of the cross section within the 99%, 95% and 90% CL. The lines are theoretical 
predictions of an enhancement of the neutrino-nucleon cross-section by electroweak sphalerons [ p4[ 
(short-dashed), p-branes (long-dashed) and string excitations |2^ (dotted). Left panel: Fit to 
the AGASA and HiRes CR data. Right panel: Fit to the AGASA, HiRes, and Pierre Auger GR 
data. 


As an illustration, we have considered three models of a rapidly increasing neutrino- 
nucleon cross section based on electroweak sphalerons |24], p-branes [36| and string exci¬ 
tations |27]. (i) Electroweak sphalerons: We have used the neutrino-nucleon cross section 
induced by electroweak instantons shown in Ref. 1241, based on the neutrino-parton cross sec¬ 


tion from Ref. |^|, the latter exploiting numerical results from Ref. |^]. A fit of the injection 
spectrum and the red-shift evolution of the source luminosity using this cross section repro¬ 
duces a 98% CL of this model, which is in very good agreement with Fig. (ii) p-branes: We 
have calculated (t{vN —> brane) given in Ref. |@ as Eqs. (9) and (10) for m = 6 extra spatial 
dimensions, a fundamental scale of gravity Mjy = 300 TeV and a ratio L/L* = 0.005 of small 
to large compactification radii. We let all partons interact universally with the neutrino. Our 
fit of 7 and n gives a 83% CL. (iii) String excitations: For the neutrino-quark cross section 


given in Ref. [27| with a string scale M* = 70 TeV and for the set of parameters Nq = C = 16, 
characterizing the width and the absolute normalization, respectively (cf. their Fig. 2), we 
derived the neutrino-nucleon cross section. Our fit gives 84% CL for this cross section, again 
in very good agreement with Fig. 


In a separate fit we have also included the preliminary CR data reported from PAO |53]. 
The result is shown in the right panel of Fig. The rate of horizontal events from PAO will 
also have a great impact on the analysis, but preliminary results are not yet published. The 
combined vertical spectrum of AGASA, HiRes, and PAO is still compatible at the 95% CL 
with our scenario as one can see in the right panel of Fig. 












































5. Conclusions 


Our analysis shows that ultra high energy cosmic rays measured at AGASA and HiRes might 
be composed of extragalactic protons and strongly interacting neutrinos. The enhancement 
of the neutrino-nucleon cross section has to be rapid in order to agree with the experimental 
results from horizontal events at AGASA and contained events at RICE. We derived allowed 
regions of a flexible parameterization of the cross section using a goodness-of-fit test and 
showed the compatibility of our results with theoretical predictions from sphalerons, p-branes 
and string excitations. The good agreement between these fits with our simple parameteri¬ 
zation indicate that out results (in particular Fig. can be used as a quick test for strongly 
interacting neutrino-nucleon cross sections. 

Our results are based on the assumption that ultra high energy protons and neutrinos 
originate at optically thin sources. One should keep in mind that a more general injection 
spectrum of cosmic rays or optically thicker sources would certainly alter the prediction of 
the neutrino fluxes and accordingly our results for the range of neutrino-nucleon interactions. 
As a rule of thumb cosmic accelerators always produce cosmic rays, neutrinos and gamma 
rays with comparable luminosities. Also the top-down production of cosmic rays is always 
accompanied with high energy fluxes of other particles. Hence, the search for the origin of 
the highest energy cosmic rays involves essentially a multi-messenger analysis. It is in the 
focus of future neutrino experiments like ANITA, IceGube or the Pierre Auger Observatory 
(horizontal events) and gamma ray observatories like H.E.S.S. and MAGIC [^] to resolve 
some of these model ambiguities in the analysis. 

Notably the Pierre Auger Observatory combining the experimental techniques of AGASA 
and HiRes as a hybrid detector will have a great impact on scenarios with strongly interacting 
neutrinos. With a better energy resolution and a much higher statistics it will certainly help 
to clarify our picture on ultra high energy cosmic rays. 
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